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We discuss the existence of breathers and lower bounds on their power, in nonlinear Schrodinger 
lattices with nonlinear hopping. Our methods extend from a simple variational approach to fixed 
point arguments, deriving lower bounds for the power which can serve as a threshold for the existence 
of breather solutions. Qualitatively, the theoretical results justify non-existence of breathers below 
the prescribed lower bounds of the power which depend on the dimension, the parameters of the 
lattice as well as of the frequency of breathers. In the case of supercritical power nonlinearities we 
investigate the interplay of these estimates with the optimal constant of the discrete interpolation 
inequality. Improvements of the general estimates, taking into account the localization of the true 
breather solutions are derived. Numerical studies in the one dimensional lattice corroborate the 
theoretical bounds and illustrate that in certain parameter regimes of physical significance, the 
estimates can serve as accurate predictors of the breather power and its dependence on the various 
system parameters. 



I. INTRODUCTION 



The discrete nonlinear Schrodinger (DNLS) model constitutes a ubiquitous example of a nonlinear dynamical lattice 
with a wide range of applications, extending from the nonlinear optics of fabricated AlGaAs waveguide arrays as in 
Q , to the atomic physics of Bose-Einstein condensates in sufficiently deep optical lattices analyzed in . Partly also 
due to these applications, the DNLS has been a focal point of numerous mathematical/computational investigations 
in its own right, a number of which has been summarized in and is related to models used in numerous other 

settings including micromechanical cantilever arrays and DNA breathing dynamics [isj . among others. 

In this work we consider a variant of the DNLS equation of the following form: 

N 

i^„ + e(AdV')« +a^nXl('^''^)»6Z" +/3|V'nP'''/'n = 0, (1.1) 

J=l 

on a A'^-dimensional lattice which can be finite if supplemented with Dirichlet boundary conditions, or infinite (n € Z^) . 
In (jl.ip . e > is a discretization parameter e ~ with h being the lattice spacing, and {Adip)n stands for the 
A'^-dimensional discrete Laplacian 

(AdV)„ez« = J2 ~ 2iV^„, (1.2) 

where Afn denotes the set of 2N nearest neighbors of the point in with label n. The nonlinear operator Tj is 
defined for every ipn, n — (rii, 712, . . . , n^v) € Z^, as 

(7j'V')«GZ" = IV'(ni,Ti2,...,rij+l,Tij + i,...,n„)P + I V'(ni .nj l,nj+i ,njv ) P ' j = ^, ■ ■ ■ (1.3) 

The nonlinearity parameters a, /3 G M. In the case a = 0, /? 7^ 0, one recovers the classical DNLS equation with power 
nonlinearity. The case where a,/? 7^ 0, corresponds to the DNLS equation with nonlinear hopping terms. The DNLS 
equation (jl.ip , is a Hamiltonian model with a Hamiltonian of the form: 

N +00 „ 

H[^]=e(-AdV,V')2-^ J2 l^(«i,«2,...,«„«,+i,...,n„)nV'(„i,„„...,„,+i,„,+,,...,„„)l'-^ ^ (1.4) 
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Let us note for convenience discuss the 1 — D lattice, where the equation reads: 

ii^n + e(V'«-l - aV'ri + tpn+l) + a^^n (|^n+lP + \ijn^l\^) + P\i^n?" i^n = 0, (1.5) 

with the Hamiltonian 

riG^ n^Z n£Z 

The Hamiltonian (|1.4p and the power (or norm) 



are the conserved quantities of this lattice dynamical system. 

We will present theoretical and numerical results related to the existence of time periodic (standing wave) solutions 
of the form 

Mt) = e'^^V™, f^eK. (1.8) 

The physical interest in this particular model stems from various contexts, as the modeling of quantum lattices 
and w aveg uide arrays and the approximation of the dynamics Klein-Gordon (KG) and Fcrmi-Pasta-Ulam (FPU) 
chains [16l - l2ll |. Eq. (11.51) for cubic (cr = 1) nonlinearity corresponds to the classical limit of the quantum DNLS 
equation introduced in [l6|. In the quantum lattice introduced therein, the inclusion of the nonlinear hopping term 
allows a fast energy propagation as long as a is high enough with respect to /3. Such terms (the additional ones 
to the classical DNLS with cubic onsite nonlinearity and linear coupling between sites) have appeared in physical 
considerations within the modeling of waveguide arrays [3, [2^ , establishing that in the case of large penetration 
length or closely spaced waveguides these terms are not negligible; however, it should be noted that in this case 
additional terms of the same (cubic) order should be included in the relevant modeling jisl [20j . Nonlinear hopping 
terms appear also from FPU and KG chains of anharmonic oscillators coupled with anharmonic inter-site potentials, 
or mixed FPU/KG chains. The generalized DNLS system of [l7j involving, among others, the nonlinear hopping terms 
considered therein has been derived as a perturbation of the integrable Ablowitz-Ladik system, by the rotating wave 
approximation on the FPU chain. A similar DNLS system has been derived in (2]| , approximating the slow dynamics 
of the fundamental harmonic in the Fourier series expansion of discrete small amplitude modulational waves. Relation 
of such DNLS systems as models for the energy transport in helical proteins has been discussed in [2^ . However, it 
is worth remarking that additional terms should also be taken into account therein, as well. Furthermore, such terms 
have been studied in their own right mathematically while considering the properties of potential traveling waves 
within a generalized class of DNLS models in [2^ . 

In this work, our main scope is to derive lower bounds for the energy of discrete breathers for the DNLS system 
(jl.ip and discuss their relevance as thresholds for their existence. In this point of view, (|l.ip seems to be of particular 
interest due to the interplay and the expected competition of the nonlinear hopping and the generalized power 
nonlinearities. Extending the arguments based on variational methods [23 - l26| and the fixed point approach of [23l to 
establish the existence of solutions (jl.Sp , we show the existence of lower bounds on the power of breathers on either 
finite or infinite lattices. The bounds depend explicitly on the dimension, and the nonlinear lattice parameters, as well 
as on the frequency of the solution. They have a simple geometric interpretation visualized in Figure (TJ elucidated in 
particular by the fixed-point approach: The energy bounds can be interpreted as the radius i?crit of the closed ball 
centered at in the energy space denoted by i?(0, i?crit). Breathers do not exist in the closed ball i?(0, i?crit), and 
a non-trivial (e.g. non-zero) breather solution being in £^ \ i?(0, i?ciit) should have energy V > i^^rit- The result is of 
physical significance related to energy thresholds (where by "energy" here we mean power, or squared norm) for 
the formation of breather solutions. In particular, it indicates that for a given set of parameters, no periodic localized 
solution can have power less than the prescribed estimates. 

It should be remarked that this result is of different nature if compared with the excitation threshold phenomenon 
of [28l. [29j for discrete breather families, possessing a positive lower bound on their energy when the lattice dimension 
N is greater than or equal to some critical dimension. In the context of DNLS systems with power nonlinearity, the 
restriction for the appearance of the excitation threshold is interpreted in terms of the nonlinearity exponent as ct > , 
[29| . In this point of view, a can be considered as critical when a = and supercritical (subcritical) when a > 
(cr < jf), and the excitation threshold exists in the case cr > It is crucial to remark that the set of parameters for 
which the excitation threshold 7?.thrcsh is apparent suggests that the energy bounds i?crit are not sharp as thresholds 
for existence/no nexistence. In particular, when -Rcrit is the value derived by the fixed-point approach, it is observed 
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that i?thresh > -Rcrit, I2J, [23 • ^01 further discussions on the excitation threshold for FPU and Klein-Gordon lattices 
we refer the interested reader to jsof . 

Section II is devoted to the derivation of the estimates by variational and energy methods employed in the case 
of finite lattices, and Section III is devoted to the fixed point approach in infinite lattices. While the methods are 
applicable for both subcritical and supercritical nonlinearities, in the latter case we investigate their interplay with the 
optimal constant of the discrete interpolation inequality of [29[ and its analytical estimation proposed in |3l| (Section 
IIIB). In Section IV we perform numerical simulations testing the lower bounds as thresholds for non-existence of 
breathers with respect to the variation of the lattice parameters, while section V briefly summarizes our conclusions. 
The previous studies proved the validity of these bounds as energy thresholds for the existence of breather solutions and 
justified that there are elements of breather families (parametrized by the lattice parameters) which tend to saturate 
the theoretical bounds in the case of large and small nonlinearity exponents. Aiming to improve this prediction for 
extended parameter regimes, we consider a refinement of the lower bounds, on account of the finite localization length 
of the true breather solutions and the expectation that the main contribution to the power comes from the central and 
adjacent sites, being the most excited. To incorporate this claim in the numerical simulations, we perform a cut-off 
procedure which considers the part of the system for the oscillators occupying a unit length around the central site 
plus the adjacent to this unit length as well. This cut-off improves the capture of the contribution of the linear part of 
the system to the power, manifested in the bounds by the first eigenvalue of the linear operator. The first eigenvalue 
estimates the contribution of the linear part from below. Contrary to the estimation of the linear part in the real 
length, its unit length approximation is not negligible since the linear mode over the latter is strongly localized. This 
is reflected in the numerical simulations performed for the case of the cubic nonlinearity. These simulations reveal 
that in the weak coupling regime the bounds are getting closer to the numerical power, and in some cases provide 
its accurate prediction. This quantitative response is observed in particular versus the nonlinear hopping parameter 
a. The good behavior of the estimates indicates that the approach presented can be promising in a study of DNLS 
systems encountered in the aforementioned applications, involving the full expansion of nonlinear hopping terms being 
however of the same order. 

We conclude the introductory section, by mentioning that although our results concern both the cases of finite 
and infinite lattices the term "breather" has been used for the standing wave solutions (|1.8p in the finite case, only 
for the sake of brevity. The important issue of the localization properties of the solutions in the transition from 
the finite to the infinite lattice is not addressed in the present work. We refer to for a detailed discussion on 
the spatial decay and stability properties of the solutions when the lattice size is varied for small-amplitudes (i.e., 
near the continuum limit), as well as, for relative localization estimates. For the convergence of solutions, defined by 
constrained variational problems in finite lattices to unimodal and even profile breather solutions (centered on a site 
or between two lattice sites) in infinite lattices, we refer the interested reader to [33| . 

a. Preliminaries. For convenience, we recall from (23 . [25| some preliminary information on various norms and 
quantities, that will be thoroughly used in what follows. 

The finite dimensional problem is formulated in the finite dimensional subspaces of the sequence spaces i^, I < p < 
oo, 

£P{Z^) = {(f>eF : (/)„=0 for |||n|||>i^}, (1.9) 
where | ||n|| | = maxi<i<jv for n = (ni, 712, . . . , njv) € ■ Note that in the case of the infinite lattice 

ll'^llg < ||0||p, 1<P<<Z<00 (1.10) 

O<e(-Ad0,0)2 < 4eiV ^ |0„|2. (1.11) 

rtGZ« 

For the finite dimensional case we have that £p(Z^) = C^^-'^'+i) , endowed with the norm 




and that the well known equivalence of norms, 

m\g<m\p<i2K+i)^^m\„ i<p<g<oo, (1.12) 

holds. 

At this point let us remark on some basic facts on the eigenvalues of the discrete Dirichlet Laplacian, since they will 
naturally appear in the estimates that will be derived in what follows and have an important role in the numerical 
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Figure 1: Simple geometric interpretation of the energy lower bounds obtained by the fixed point argument: Breathers do 
not exist in the darker (red) area, the closed ball _B(0,i?crit) of , centered at and of radius i?crit. The lighter (green) area 
represents the area of the energy space where breather solutions exist. Although the non-existence result does not depend on 
the dimension and the lattice parameters, the radius -Rcrit of the closed ball _B(0,-Rcrit) of non-existence, quantitatively is a 
function of the lattice parameters a, j3,a, the frequency and the dimension of the lattice N . Note that i?crit is not sharp with 
respect to non-existence. This is suggested from the case for which the excitation threshold i?thrcsh is present. In this case it 
is possible that iithresh > Rciix , and the dark (red) area is enlarged. 



simulations. For the lZ?-latticc of K + 2 oscillators, ri = 0, . . . , i^T + 1, let us consider the discrete eigenvalue problem 

for (f)n e M, 

- eAd(t>n = ^J■4'n, n = l,...K, (1-13) 
with Dirichlet boundary conditions, (f)o — (f'K+i ~ 0. Starting from the standard case 

-TT where h = — — ^ . 
where L denotes the length of the chain, the eigenvalues are 



e = — r where h=— — -, (1-14) 



^„(/i) = — sin^ — ^^^^sin2 — -— , n = l,...,A'. 



\2L J \2{K + 1 

Thus, in the case (|1.14p the principal eigenvalue is 

riiih) = ^sm"^ — = ^-^sm^ — . (1-15) 

' \2Lj \2{K+1)J ^ ' 

The discrete system is modeled when h ~ and in the limits /i — > and /i — > oo we have 

lim = Ai = — -, (continuous limit), (1-16) 

lim /ii(/i) = 0, (anticontinuous limit). (1-17) 

h— f CXD 



In the particular case of L = 1 we have 

lim Ui(h) = / 

4 < mih) < TT^ for 0<h<l. (1.19) 



lim fii{h) = Ai = TT^, (1.18) 
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In a general discrete case the parameter e > can be either related or not related with the lattice spacing h. As an 
example for the former, we may fix the linear coupling constant e > 0, varying the number of oscillators, equidistanced 
with lattice spacing h = j^^- We have 

= 4esin2 f ^ ) = 4esin2 f — - ) , lim = 0, {h ^ when K ^ oo) (1.20) 

\2L J \2[K + Ijy /i— i-o 

0<Mi(/i)<4e. (1.21) 

Increasing K, (ll.20p - p.21l) can be considered as a particular approximation of an infinite lattice. Note that in the 
case of the infinite lattice Z^, for the discrete Laplacian with e = 1, wc have that cr(— A^) C [0,4A^]. 

Relations ()1.15|) and ()1.16|) . (|1.17|) arc valid for a general coupling (depending or not depending on the lattice 
spacing) behaving as e ^ with e sufficiently large. Similar observations arc valid in the case of the N-dimensional 
discrete Laplacian. 

Finally, we recall that the variational characterization of the eigenvalues of the discrete Laplacian in the finite 
dimensional subspaces £'^{'E^), showing that ni > 0, can be characterized as 

Ml = mf =^ rrr?- (1-22) 



^0 



Then, (|1.22p implies the inequality 



/ii l-^nl' < e(-Ad,^,</))2 <4eiV ^ |0„|2. (1.23) 

lll"lll<-ff lll«lll<^ 



II. FINITE DIMENSIONAL LATTICES 



This section is devoted to the DNLS equation with nonlinear hopping terms a,/3 7^ 0, supplemented with Dirichlet 
boundary conditions 

N 

ii^n + e(AdV)n + XlC^'^)"^^;" + /3|V'nP''V'n = 0, (2.1) 

V-n = 0, ||n|| > K. (2.2) 

We will employ a constrained variational approach on the nonlinear energy functional involving the nonlinear hopping 
term. Noticing that the existence result can be established by minimization of the Hamiltonian or by application min- 
max methods (e.g mountain pass type theorems), the usage of alternative functionals may reveal interesting conditions 
on the nonlinearity parameters. An example is given in |26l . Section 2.2 & 2.3, pg. 9-18], where the minimization of a 
linear energy functional under a nonlinear constraint verified conditions for the co-existence of breather profiles. For 
instance, this alternative approach for (j2.ip will show the existence of a regime for the hopping parameter a where 
an upper bound for the power is valid (see Remark lII.3|) . 

Note that the case of Dirichlet boundary conditions is of interest in particular for numerical simulations; since the 
infinite lattice cannot be modeled numerically, numerical investigations should consider finite lattices with Dirichlet 
or periodic boundary conditions. The latter should be imposed for moving breathers colliding with the boundary. We 
expect that the variational approach can be applied in the case of periodic boundary conditions, but the details have 
to be checked. 

We shall consider first the focusing case for the parameters a, /S > and we shall briefly comment on the dcfocusing 
one a,P < which can be treated similarly. 

A. The focusing case a, 13 > 0-Solutions ^„(t) = e'"Vn, > 
Substitution of the solution (jl.Sp into (|l.ll) shows that 0„ satisfies the system of algebraic equations 

N 

- e(Arf0)„ + r!(/)„ - a(/)„^(7^-0)„gz« -/Sl^nP'^n = 0, e R, \\n\\<K, (2.3) 
0„ = 0, \\n\\>K. (2.4) 
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Let us note that in the anticontinuous Umit e = 0, the corresponding energy equation reads as 

N 

\\n\\<K \\n\\<K j=l \\n\\<K 

Its positive right-hand side, imphes direetly that in the hmit e = 0, the foeusing case supports only solutions with 

n>o. 

For e > we will also restrict our considerations to the case of solutions with i7 > 0. We recall two auxiliary 
lemmas regarding the differentiability of the nonlinear terms if viewed as nonlinear functionals, which can be proved 
as in [23, Lemma 2.3, pg. 121]. 

Lemma II. 1 Let <f> E . Then the functional 

N +00 

ni ,n2,...,nj ,nj + i ,...,njv) I I */'(ni ,712 ,. . . ,nj +l,nj + i ,. . .,n7v) I ' 

j — 1 nj — — 00 

is a C^(£^,R) functional and for all ijj £ t'^ , 

N +00 

< "^''(0)iV' >= 2Re^ ^ l<?^(ni,n2,...,ni+l,ni+l,...,nN)P'^("l,n2,...,n3,n3 + l,...,"N)V'(,ii,n2,...,nj,nj + i,...,riN) 
j — 1 rLj— — oo 

N +00 

-|-2Re ^ ^ I0(ni,n2,...,nj-l,nj + i,. ..,««) I ,"2 ."j + l , ■ ■ ■,nN)''/'(Tii ,n2 ,. . . .rij ,nj + i , . . .,rijv) • (2-5) 
j — l nj— — oo 

Lemma II. 2 Let G Then the functional 
is a C^(£^,M) functional and 

</:'(0),V>=2(a + l)Re ^ |0„p-(/)„V^„. (2.6) 

Both Lemmas IIL II and IIL 2 1 remain valid in the case of the finite lattice (space £^(Z^)). 

The first result on the existence of time-periodic solutions (|1.8p of (|2.ip . is via a constrained minimization problem 
for the functional 

:=e(-Ad(/),0)2+r! ^ |0„p-aV(0), Q > 0, a>0. (2.7) 



Theorem II. 3 A. Consider the variational problem on ^ (Z^) 



inf : =m| 



(2.8) 



for some 57 > 0. Then, there exists a minimizer (j> G £^(Z^) for the variational problem \2. and I3{M) > 0, both 
satisfying the Euler- Lagrange equation i2.S\) - [K^ and X^nGZ™ l^nP'^"''^ = M{a + 1). 

B. Assume that the power of a solution of the problem i2.3\) - (K^ is Vlcj)] — . Then the power satisfies the lower 
bound 

Rlf<R^^V[$], (2.9) 
where R*j denotes the unique positive root of the algebraic equation 

/3x'" + 2aiVx' - (Ail + f7) = 0. (2.10) 
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C. We assume that 



a > 1 

Then a breather solution of lll.l] } satisfies the lower bound 



1 

2P 



f2 + /ii 



in either one of the cases 

(i) (lattice spacing condition) For all VL > Q if 



(ii) (frequency condition) For all e > if 



e > 



n > 



{2aN)^ cr - 1 
(/3a) ^ ^ ^ 

{2aN)^ cr - 1 
(/3a) ^ -^ic^ ' 

{2aN)^ a - 1 



(2.11) 



(2.12) 



'— a 



(2.13) 



(2.14) 



Proof: A. Wc consider the set 



1 



a + 1 



From Lemma [2T5l we may easily infer that £ : ^ R is a C^-functional. Moreover, by using inequahty (jl.l2p . we 
deduce that 



> -aN 



l>~aNU\\i 



> -aN{2K +1)^{£[(I)])^ 

= -aN{2K + l)^{M{a + l))^. 

Therefore, the functional £ : — >■ M is bounded from below. By the definition of the set B„ and the fact that we 
are restricted to the finite dimensional space ^^(Z^), it immediately follows that any minimizing sequence associated 
with the variational problem (j2.8p is precompact. Hence, by the Weierstrafi minimization theorem [341 Proposition 8, 
pg. 37], any minimizing sequence has a subsequence converging to a minimizcr and £ attains its infimum at a point 
(f) in Ba- To derive the variational equation (|2.3p . we consider first the C^-functional (due to Lemma |TL2]) 

1 



and we observe that for any (p € B^ 



CM[^] = —-m-M, 

a + 1 



{C',j[^],cl,)=2m^2M. 



Thus, the regular value Theorem ([H, Section 2.9], [H, Appendix A,pg. 556]) implies that the set B^r = £^/(0) 
is a C^-submanifold of £'^{'E^). Application of the Lagrange multiplier rule, implies the existence of a parameter 

l3 = l3{M) e R, such that 



£'[^]-f^jr'Mi; 



2e(-Ad(^,i^)2 + 2f]Re ^ (^„V„ 

AT +00 

2aRe ^ ^ ^ ] I0(ni,n2,---,n3+l,nj-|-l ,■■■,"«) I 0(ni ,Tt2 ,■ ■ -."j ."j + l ,■■■."«) '/'(ni ,n2 ,.. ..n^ + i ,...,««) 
j — 1 rij — — oo 

N +00 

2aRc y ^ y ^ I0(ni,n2,...,ii3-l,n3 + i,...,njv) I '^(ni ,n2 ,.. .,nj .Mj + i ,...,««) ,712 -."j ."j + l ,■■ 

j — 1 n J — — 00 

(2.15) 



2/3Re 



,p-0„^„ = 0, for aU V e ^^{^k)- 
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Setting -0 = in (|2.15|) . wc find that 

nGZ" 

JV +00 

— 2aRe ^ ^ l0(ni,n2,---,nj+l,nj+l,---,raN)l I ,"2 .ij + l ,■ ■ ■,"«) I 

TV +00 

— 2aRc ^ ^ |0(ni,n2,...,nj — l,nj + i,...,njv) I I '^("l ,"2 ,■ ■ -."j ."j+l ,■■■,"«) I 

j^l rjj— — oo 

= /? E l^nl^'^''^. (2.16) 

By virtue of (|1.23p . we deduce that the following estimate 

riGZ" 

> ^ilml + ilM\l-2aNmi (2.17) 
holds. Let us assume that Vlcj)] = \ \4>\\2 = ■ Then from p.l7p . we obtain that 

F[<j)] > R^ipi +n- 2aNR^). 



Therefore, assuming that 



or assuming in terms of a that 



we deduce that > 0. Since cj) € B^j cannot be identically zero and JF[(I)] > 0, it follows from p.l6p that /3 > 0. 
Summarizing, we have proved that for given J7 > 0, there exists a minimizer (j) and a Lagrange multiplier /3 > solving 
the variational equation p.l5|) . Clearly a solution of the variational equation (|2.15p is a solution of the Euler-Lagrange 
equation 

B. It is necessary to verify first that any solution of (|2.3p - (|2.4p is a solution of the minimization problem (|2.8p . 
Indeed, if is a solution of p.3p - (|2.4p . multiplying p.Sp by in the P{'L^) and using the Dirichlet boundary conditions 
we infer that satisfies equation (|2.16p . written as 

F[4>] = PC[4>]. (2.20) 

Then, due to Lemmas III.ll and 111.21 solves also the equation 

(j-'[0],V) =/3(/:'[0],V), for all i^efiK)- 
Comparing (j2.15p with ()2.16p it can be easily seen that the equation above is equivalent to 

(^E'[4>1^)=P{C'[4>1^), for aU V> e (2.21) 

thus, is a minimizer of the minimization problem (j2.8l) . The converse follows immediately by (j2.2ip and the fact 
that in the discrete setting a "weak solution" of (j2.2ip coincides with a solution of (|2.3p - p.4l) . Furthermore, by setting 
i/; = in (|2.2ip we recover that satisfies the equation (|2.20p . 

Assuming now that the power of the solution of (|2.3p is V[(j)\ ~ H^Hl ~ i by using (|1.12p and (|1.23p we get from 
(|2.20p . that R satisfies the inequality 

fii+fl<2aNR'^+PR^''. (2.22) 
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The algebraic equation p.lOp considered for x S [0, oo), has exactly one positive root < i?*./. Then, comparison of 
the equation (|2.10p with inequality (|2.22p . implies that the power V[<t)\ must satisfy the lower bound (|2.9p . 
C. Applying Young's inequality 

ab < -aP + -^-j-b"^, a,b> for any e > 0, 1/p + l/q = 1, 
p qeiiP 

with p = (7, q = a = , b = 2aN and e = /3cr we get that 

2aE? < BR^" + 4 (2.23 

{Pa)— ^ 

Inserting (|2.23p into (|2.22p we derive the lower bound (|2.12|) . o 

Remark II. 3 1. (The lower bound for the cubic nonlinearity) For the case of cubic nonlincarity cr = 1, inequality 
(|2.22p implies that the power of the periodic solution V'ri(i) = e'^^'c^n, > must satisfy the lower bound 



2aiV + /3 



<R^^V{4)\. (2.24) 



2. (Interpretation of condition ()2.18p ). The result of Theorem (jll.3| establishes for arbitrary given f2 > and a > 0, 
the existence of a nontrivial (/> e i!^(Z^) and the existence of /3 > as a Lagrange mulitplier such that V'n(0 — e'^^*(/>n, 
solves equation ()2.ip with /3 > as a parameter for the power nonlinearity. On the account of this result, the meaning 
of condition (|2.18p is that there exists /3 > and a range of the hopping parameter < a < a* for which the associated 
minimizer (j) has power satisfying the upper bound 



2 ^ Ml + ^ 



Note that the existence of the range of the hopping parameter a stated above is also established by (|2.18p -see (|2.19p . 
3. (Case a — ?> 0, /3 > 0-DNLS with power nonlinearity). The proof of Theorem (jll.3p remains valid for the case a = 0, 
where one has to consider the constrained minimization problem p.Sp for the functional setting a = 0. Thus for 
the classical DNLS with power nonlinearity we recover from inequality p.22p . the lower bound 



Ml + 



<R'^V[<j)]. (2.26) 



The lower bound (|2.26p is the same as (5.27) and (5.31) of 0] for the DNLS with power nonlinearity. 

B. The defocusing case a, /3 < 0-Solutions Vn(t) = e"'"Vn, > 0. 

We shall briefly comment on the existence of breather solutions, for the case of negative nonlinear parameters 
q;,/3 < 0. We set for convenience a = —k, f3 = —X where k,X > 0. It should be remarked that the case of negative 
parameters can be reduced to the case of positive ones, under the staggering transformation. We recall that this 
transformation is defined as 

N 

^„ -> (-l)l"l^„, |n|=^7i,, (2.27) 

i=l 

(see e.g. the discussion of [s^, pg. 7]). The case of negative parameters, corresponds to the existence problem for 
solutions 

^n{t) = e^'^Vn, > 0, (2.28) 
where 0„ satisfies the system of algebraic equations 

N 

- e(Arf(/.)„ - r!0„ + K(/)„^(7^-(/.)„ez" +A|0„P"0„ = 0, H > 0, \\n\\<K, (2.29) 

j=i 

0„ = 0, ||n|| > K. (2.30) 
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The proof of the existence of breather solutions (|2.28|) is very similar to that of Theorem IIL3[ and we refrain from 
giving the details. We just note that the constrained minimization problem will consider the C^-functional 



J2 + kV(0), n>o, K>o. 



(2.31) 



Theorem II. 4 A. Consider the variational problem on £'^{Z^) 

infirm : =m| 

for some O > 0. Assume further that 

n > ieN. 



(2.32) 



(2.33) 



Then, there exists a minimizer (j>* G £^(Z^) for the variational problem S2. 32]) and X{M) > 0, satisfying both the 
Euler- Lagrange equation ^MB-(KW o.nd E„gz" ICP"^^ = 

B. Assume that \2. 3S\] holds and that the power of a solution of the problem i2.29\) - ^2.S0\) is Vl 
power satisfies the lower bound 



where R^^d denotes the unique positive root of the equation 

XX^" + 2kNx^ - (rj - AeN) = 0. 

C. Let cr > 1 and assume that 

Then the power satisfies the lower bound 

n - ieN - 




*] = Then the 
(2.34) 

(2.35) 
(2.36) 



(2.37) 



Remark II. 5 1. (The lower bound for the cubic nonlinearity) For the case of negative parameters a 



-A, 



At, A > and of cubic nonlinearity (7 = 1, the power of the periodic solution ijjnit) = e (j)m ^ > must satisfy the 
lower bound 



2kN + \ 



< R^ = V[(t)*], n > AeN 



(2.38) 



2. (An upper bound for some range of parameters). The result of Theorem III. 41 establishes for given > 4e and 



a = —K < 0, the existence of a nontrivial (f)* € £ (Z^) and the existence of /3 = — A < such that ijjn{t) 
solves equation p.ip with /? > as a parameter for the power nonlinearity. As in remark [II. 3I -2. a similar condition 
to (j2.18p can be derived, implying that there exists a parameter A and a range for the hopping parameter k for which 
the corresponding minimizer cj)* has power satisfying the upper bound 



O _ AfM 



(2.39) 



3. (Case k — > 0, A > 0-DNLS with defocusing power nonlinearity). The proof of Theorem III. 41 remains valid for the 
case K = 0, where one has to consider the constrained minimization problem (|2.32p for the functional £, setting k = 0. 
Thus for the classical DNLS with power nonlinearity we recover the lower bound 



n - AeN' 



A 



<R^ =V\ 



n > 4eN. 



(2.40) 



The lower bound (|2.26p is exactly the same with that derived in |2J] for the one dimensional DNLS with defocusing 
power nonlinearity. 

4. Condition (j2.33p is related with the extension of the phonon band for defocusing-type DNLS equations, to the 
interval [0, 4eA^]. Combining the results of Therorcm III . 31 for the focusing case and of Theorem III. 41 for the defocusing 



one, wc have that for breathers in the ansatz tp^i 
(defocusing case) and 17 < (focusing case). 



J, frequencies f2 G R, must lie in the intervals f2 > 4eA^ 



11 



III. INFINITE Z^, N >1 LATTICES 

For the infinite lattice Z^, we will consider the problem of energy bounds for breathers of the DNLS (jl.ip by a 
fixed-point method. The method establishes that the stationary problem (|2.3|) defines a locally Lipschitz map on the 
phase space When the map is a contraction, gives rise only to the trivial solution. The Lipschitz constant for 
the contraction mapping defines the critical power above which we should expect existence of breathers. Below this 
critical power there is non-existence of breather solutions. The Lipschitz constant contains all the lattice parameters, 
including the dimension of the lattice and the frequency of the solution. 



A. The case a, /3 > 0-Solutions i/'n (i) = e'*^*, SI > 0: Fixed point method 

The infinite system of algebraic equations (|2.3p for breathers in the case of the infinite lattice will be treated by a 
fixed point argument. We recall that the linear and continuous operator 

- eAd + n : f ^ e^, (3.1) 

satisfies the assumptions of Lax-Milgram Theorem [H, Theorem 18. E, pg. 68], since 

e{~Ad<l>,<l>)2 + nm\l>nm\l for all ^ef. 

This is the first step to verify that for given 2 € the auxiliary problem defined by the linear operator equation 

N 

- eArf(^„ + n(j)n = aZn^{TjZ)neZ« + P\Zn\^'' Zn, (3.2) 

i=i 

has a unique solution (j) E P . The second step, according the Lax-Milgram Theorem is to justify that the right hand 
side of (|3.2p is in 1°^ \i z ^ P . Indeed, by using the inequality 

p 

^ I'/'nT < ( ^ , foraU l<g<p<oo, (3.3) 

for p = 4(7 + 2 and q = 2, it follows that 

Furthermore, for the nonlinear map J : IP' P' ^ 

N 

J'[Zn] = y~^(7^-^)ngZ"7 

we have 

\\J[z]\\l < 2N sup |z„|2 V \z„\^ < 2N\\z\\i 

Therefore we are allowed to define the map A : i'^ ^ P , by A{z) := 0, where is a unique solution of the operator 
equation p.2p . Clearly the map A is well defined. Let C,^ be in the closed ball 

BR:={z(^e :||z||p<i?}, 

and (f) ~ ^(C), "0 ~ -^(0- The difference x := — 0! satisfies the equation 

- eAdXn + ^a{J[Cn] ~ J[Cn]) + Pi\Cn\^'' (n " \^n\^'' ^n) ■ (3.5) 

We consider the linear and continuous operator M : ^ 

N 

[Zn\ ^ ^ ^(ni ,712 ,...,nj +l,nj-|-i , . . ..n^v) ~t~ -^(ni ,n2 , . . — l,nj-f i ,...,n7v) i j ^ ; • ■ • : ^ 
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satisfying 

\\M[(j)]-MM\2<2N\\(j)-^p\\, for all € (3.6) 
Then, the first term of the right-hand side of p.Sp ean be written as 

By using p.6p and inequality p. 31) for p = 4 and q ^ 2, we observe that 

< sup^|A^[|C„|^]|^ I?»-C«l' 

2 



< 47V^('E E 1^" 



< 47V^i?^||^-Cll2- (3.7) 

Using again (|3.7p we get that 

||C(A<[|^lV-^[ICn)ll^ - E \Cn\'\Mmn\']~M[\Cn\T 

n6Z» 

< AN' sup |C„P E ll^-l'-ICnH 

< 4iV2i?2 sup (|e„| + |C„|)' E l?"-Cn|' 

< 8iV2i?4||^_^||2. (3.8) 

Hence, from (|3.7p and (|3.8p . the inequality 

||^K]-^[C]||2<%/l2A^i?'||e-CI|2 (3.9) 
readily follows. Moreover, it holds that (cf. [U Lemma II. 2]) 

E llCnP^Cn - l^nP^CnP < (2a + 1)' R^'^ E 1^" " ^"l"' (3.10) 

nGZ" nSZ" 

Now, taking the scalar product of p.Sp with x in and using (|3.9p and (|3.10p . we have 

e(-A,x, x)2 + f^llxll^ < c,\\xMJm - J[C]\\2 + PWxM ICI^^C - le'-Ch 

< LiR)\\xMC-^\\2, (3.11) 

where 

L{R) = VUaNR^ + P{2a + 1)R'''. 
Since (— A^x, x)2 > 0, from p. lip we get the inequality 

m\\l<^\\C-m + ^\\x\\l (3.12) 

From p.l2p . we conclude that 



^.\\A{z)-Am\i<^^\\c-^\\2, 



2 
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and, hence, the map A : -B_r — > -B_r is Lipschitz continuous with the Lipschitz constant 

M,«) . 

The map ^ is a contraction, and hence, has a unique fixed point if 

M{R) < 1. (3.13) 

This unique fixed point is the trivial one, since .4(0) = 0. Wc consider the polynomial function 

Il{R) L{R) ~ n. (3.14) 

The threshold value for the existence of nontrivial breather solutions can be derived from condition (|3.13p . as in the 
proof of Theorem III. 3B : Denote by i?crit the positive root of the polynomial equation n(i?) = 0. Then n(i?) < for 
every R e (0,i?ciit), that is, condition p.l3p is satisfied ii R £ (0, i?crit)- Therefore breathers of arbitrary energy do 
not exist. A breather should have power R^ > i?^rif summarize in 

Theorem III.l We assume that the parameters a,/3,a > Let Rent > denote the unique positive root of the 
polynomial equation Il{R) = 0, where n(i?) is given by |j'.-?^[ ) Then a breather solution ipn{t) = e'^*(/>„, for any > 
of mW must have power V > i?^^;^. 

The simple geometric interpretation of Theorem IIII . 1 1 is visualized in Figure [TJ Breathers do not exist in the sphere 
B(0,i?crit) of the energy space 



B. Estimates for supercritical nonlinearity exponents a > 2/N. 

A different version of dimension-dependent estimates in the case of the infinite lattice can be produced by using 
the discrete interpolation inequality of [29| 

E l'/'"l"^+'<^* f E l-^"!') (-Ad0,<^)2, a>|. (3.15) 

However, since (j3.15p is valid only for a > N/2, the derived estimates will refer only to this range of parameters. We 
recall that the range <j > N/2 is related to the appearance of the excitation threshold for breathers on DNLS lattices 
with power law nonlinearity. 

We start by multiplying (|2.3p by (p and summing over Z^, to get the equation 

N 

e(-A,</),</))2+f] E l'/'«P=a E El^»l'('^-^)»ez«+/3 E IM'^^'- (3-16) 
Using p.lSp in order to estimate the {—Ad4>, 0)2 term of (|3.16p we have 

E ^ " E Ei'^nP(7^-'^)n.z«+/? E 



< 2aN(j2 1-^"!') E 1-^" 



2CT+2 



(3.17) 



The inequality p.l7p can be rewritten as 

nR' < 2aNR* +(^(3- ^-^^^ I'^"!'"'''- (3.18) 

By using p.3p . this time for p = 2cr + 2 and q — 2, the term X^nez™ of (|3.18p can be estimated in terms of 

the power J2nez« = -R^, as 



E \M'^^'<( E i^«n 



32cr+2 
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Thus, from p.lSp and the above estimate, we derive that 



nR^ < 2aNR* + ^ - 



R 



2cr + 2 



implying that the power satisfies the inequality 



n + <2aNR^ + fiR^" . 



(3.19) 



Theorem III. 2 Assume that a > 2/N and the parameters a,/?, 51 > Let Rent > denote the unique positive root 
of the polynomial equation 



2aNR^ + PR^ 



c. 



0. 



Then a breather solution ipn{t) = g 4>n, for any ft > of must have power V > R, 



2 

crit ■ 



For an even more explicit estimate, at least an estimation of the optimal constant C* is needed. This is provided 



by 



Proposition III. 3 Let a > 2/N. There exists i^crit > 1/2 such that the optimal constant of the inequality 113.15]) 
satisfies 



^<C.<-^^^^, iV>l. 
47V ^2iy,,it - 1 AN 



(3.20) 



Proof: One of the fundamental results of ,29| is the characterization of the optimal constant C* involving the 
excitation threshold for breathers of the focusing DNLS equation with power nonlinearity. For instance it is known 
that 



thresh 



(f7 + l)e 



On the other hand, it was proved in [sTl, Proposition II. 1, pg. 6], that there exists z^crit > 1/2 such that 

< i?thrcsh < [4e7V(cr + l)]^ . 



V^^crit - 1 4iVe(a+l) 



t'crit 



2ct + 1 



(3.21) 



The estimate p.20p follows by inserting the characterization for 7?.thrcsh into p.2ip . o 
Together with Proposition IIII.3I Theorem IIII.2I can be restated and refined as follows. 

Theorem III. 4 We assume that 

a >2 when N = I and a > 1 when N > 2. 
Then a breather solution of satisfies the lower bound 



1 



AeN yJ2v„;t - 1 {2aN)— a - 1 



2(7+1 Ve 



(/3a) 



< R^ 



(3.22) 



(3.23) 



in either the cases 

(i) (lattice spacing condition) For all Q, > if 



e > 



(2aAr)^ (cr- l)(2cr + l) i/crit 



{13a) — 



ANa 



V2l'crit - 1 



(a) (frequency condition) For all e > if 



n > 



{2aN)^ cr- 1 
(/3ct)^ cr 



(3.24) 



(3.25) 
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Proof:. Inequality p.l9p can be strengthened from below by replacing 1/C* by its lower estimate as indicated from 
(12201) • Then, ([5?^ comes out exactly as in Theorem HOl C. o 

We remark that in the case of the limit a = 0, if we will repeat the calculations leading to the energy equation 
(|3.16p and inequalities p.l7p - p.l8p . we derive the inequality 



< ^IR^ < (3 



(3.26) 



Now. the positivity of the right-hand-side of p.26p implies that in the limit a = 0, the f7-indepcndent lower bound 



< R^ 



is satisfied. 

Let us also remark that the non-existence result of Theorem IIII.ll is valid in finite lattices, due to the validity of 
inequality in the subspacc oil'^{1^). Thus, the resuh can be proved in the case of finite lattices without 

any additional implications. Similarly, inequality p.l5|) is also valid in and the estimates of Theorem IIII.4I 

can be proved to be valid in finite lattices. The estimates of Theorem IIII.4I for the case a > 2/N, will be tested 
numerically in the next section. 



IV. NUMERICAL STUDY 



We present in this section, numerical results testing the behavior and relevance of the theoretical estimates, in the 
case of the ID lattice. The structure of this section has as follows. In Sec. IIV A ll we analyze theoretically a refinement 
of the original variational estimates on the example of the focusing case a, /3 > 0, aiming to improve the capture of the 
contribution of the linear part of the system to the power. This contribution is manifested in the bounds, by the first 
eigenvalue of the discrete Laplacian. The refinement takes into account the localization of true breather solutions, by 
performing a "cut-off" procedure, focusing on the most excited states. The improvement is reflected in the numerical 
simulations performed in Sec. IIV A 21 for the case of the cubic nonlinearity a = 1, showing in particular, that in some 
cases of the weak coupling regime, the estimates provide an accurate prediction of the numerical power. In Sec. IIV A 31 
we present the numerical results for the case of the quintic nonlinearity a = 2. The refined variational estimates are 
valid, due to the translational invariance of the "cut-off" procedure, even in the case of the infinite lattices, and have 
been tested against the interpolation estimates (e.g. those by the interpolation inequality of Gagiiardo-Nirenberg 
type). It was interesting to observe that the refined variational bounds give a better qualitative prediction when the 
nonlinearity parameter (3 is varied, while the interpolation estimates behave better for large values of frequencies Q. 
Finally, in Sec. IIVBI we present an indicative numerical study of the interpolation estimates in the defocusing case 
a < 0, /3 < 0. The main finding here is that the theoretical predictions are improved for large values of the parameters 
and a. 

We note that in all the numerical simulations, the results have been obtained for a ID-lattice of K = 101 particles. 

A. Focusing case (a > 0, /3 > 0) with Dirichlet boundary conditions. Solutions il>n{t) = e'^'(/)„, SI > 0. 

1. Theoretical analysis of the "cut off" procedure. 

According to the results of Theorem III. 31 A., without any restrictions on the exponent cr > of the nonlinearity, 
the first lower bound comes from the positive root i?* ./ of the equation (|2.10p 

Px^" + 2aNx^ -{fii + n) = 0, CT > 0, N>1 (4.27) 

Then any breather solution has power Vlcj)] satisfying the lower bound 

Rlj < V, for all cr > 0, N>1 (4.28) 

In the particular case of the cubic nonlinearity this lower bound reads as 



16 



Due to its relevance from a physical point of view, we have chosen the cubic nonlinearity for a first numerical test. 
The principal eigenvalue in (j4.29p manifests the contribution of the linear part of (jl.Sp . The variational characterization 
of the principal eigenvalue (|1.22p , shows that the contribution of the linear part to the real breather is estimated from 
below by the eigenvector (j)^ corresponding to the principal eigenvalue /^i, since the infimum in (jl.22p is attained by 
6^ as 



(-eArf0\0i)2 



(4.30) 



lll"lll</^ 



and (|1.22|) holds for all (j) G €^(Z;^). Qualitatively and geometrically, this approximation of the linear part seems 
reasonable, especially for breather solutions without sign changes (zero-crossings), since the eigenvector (j)^ has no 
sign-changes. On the other hand, real simulations should consider a sufficiently large chain length L, especially when 
the infinite chain is modeled in order to avoid the influence of boundary conditions. In this case, /xi — >■ (see ()1.15|) - 
(|1.20[) ') and the contribution of this approximation becomes negligible. This can be explained physically, taking into 
account the fact that the real breather solution has a localization length iioc « L while the eigenvector is extended 
through the entire chain length L. Proceeding further, since the contribution to the power outside the breather width 
iioc is also negligible, we could "cut-off" the estimation procedure, estimating the power in Lioc and the contribution 
of the linear part by the principal eigenvalue ^ixiaa of (|l-13p considered on L\oc- Practically, since the breather width 
Lioc is unknown, we may perform this "cut-off" procedure in an interval close to the interval of unit length L = \, 
expecting that the main contribution to the power comes from the excited sites included in the unit interval. This is 
certainly true for breathers centered around the center of the interval [— L, L] located at the site n = -^^j-^- It should 
be remarked that a breather can be always centered around the principal site, especially in the infinite lattice due to 
the integer translation invariance therein. 



For instance, we will consider the interval U = [— |, ^] together with the first neighbors adjacent to the points — ^ 



and i . We assume that the breather configuration is described by the vector (f) € M^+^ 



where </>„ := (/)(x„), x„ = -f 
of unit length U = [— |, |] is 



9!) = ((/)o, </)!,... (/)K+i), 00 = = 0, (4-31) 
+ nh, n = 0, . . . , K + 1. The number of oscillators located outside the piece of the chain 



9^2 



1){K + 1) 



(4.32) 



where \x~\ = min{n G Z | n > x}, a; € R. Then the number of oscillators included in the unit interval U is 

m = K + 2-6. (4.33) 

We also assume that the neighbors adjacent to the endpoints of U, are located at the sites k and fc + rn + 1. Note 
that these neighbors coincide with the end-points of U only when |; £ N. The distance y > of these neighbors from 
the endpoints of U is given by 



1 - (m - l)/i , 1 , , 
2/ = 0, if e N. 



(4.34) 
(4.35) 



We denote by U' the interval occupied by the m oscillators in U and the two neighbors adjacent to the endpoints of 
U, i.e, containing m + 2 oscillators. The length of U' is 

L' = 1 + 2y. (4.36) 

We have the following 

Proposition IV. 1 Let e>0, N = a = I. Then the power of the m oscillators included in the interval U 
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satisfies the estimate 

<'Pu<V, (4.37) 

2a + p 

where the number of points m in U is given by 

Proof: The breather configuration vector cj) in (|4.3ip can be decomposed as = ^^^'^ + (p^ 

cj)^\^ = (00, 01, ■ • ■ 0fc-l, 0fc, 0, . . . ,0, . . . , 0K+2) , (4.38) 

0^ ^ (O,...,O,0fc+i,0fe+2,...,0fc+™,O...,O). (4.39) 
Since the decomposition is hnear, at first glance the elements 4>^\^ and (fP satisfy the equations 

~e^dcl>'i + n€~o^€ (|0„+iP + |0„-iH+/3|0„|'>,^ = O, n = 0,...,i^ + 2, 

However, on the account of (j4.39p . the equation for (pu can be written as 

-e^dC^Z + f^0,V-«0,^ (10^+1 |2 + |</,^_,|2)+;3|0^|2-^^ = O, n = A; + l,...,fc + m + l, 



''fc — Vk+m+l 



= 0. 



Relabeling for convenience, the system for (ff can be considered on the interval U' of the m + 2 oscillators j = 

0. . . . , m + 2 as 

- eA,07 + 1707 - "-^7 + iCiH + /^l^yi'"'^^ = 0. J = 1, . . . ,m, (4.40) 

^li" = 0™+i = 0. (4.41) 
We also consider the linear eigenvalue problem on U' 

-e/\d4>j = t^4>j, j = l,...,m, (4.42) 
00 = 0™+i = 0. (4.43) 

The principal eigenvalue ^ii^u' of l|4.42p - (|4.43p is given by 

III U' = 4esin^ ( ) = 4esin^ ( r I . (4.44) 

^''"^ \2L'{m+l)) \2{m+l)) ^ ' 

Repeating the calculations of the proof of Theorem III. 31 on the system (|4.40l) - (|4.41l) . we derive that 

iii jji + 12 
o\r <'Pu<V, a = l, N>1, 
2a + p 

1. e, the left-hand side of (|437l) . The left hand side follows from the fact that V = YlfJl^ |0jp > Py. o 

Remark IV. 2 The estimate (|4.37p will be useful for the numerical simulations since it is valid for any e and can be 
used for the fully discrete case, even in the case of an infinite lattice, since the interval U' where the procedure takes 
place, is the same independently of the length of the chain. Thus even in the case h > 0.5, where the unit interval 
U contains only the centered site, we may perform the "cut-off" procedure for the centered site and the two adjacent 
neighbors occupying the interval U' of length L' = l + 2y = 1 + 2{h — i). The estimate (|4.37p reads as 

2alp <^^<^' ^'-''^ 
estimating the power of the breather in terms of the "most excited site" . 

In the case we approximate the continuous limit by considering e > sufficiently large, we have 
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Proposition IV.3 Let e = N = a = 1. A ssunie that e is sufficiently large, or ^ G N. The power of the 
oscillators included in the interval U 

n=k+l 

satisfies 

4(m+l)^sin^(^)+ » 
where the number of points m in U is given by {4^.33^ . 



<Vu <V, (4.46) 



Proof: Working as in tlie proof of Proposition llV. 1 1 we estimate tlie linear part of (|4.40p - (|4.4ip . by using the principal 
eigenvalue of the linear problem (|4.42[) - (|4.43|) . where in the case e = -p-, is 

t^i,U' = 7777 = 4 , , ' sm^ — — — , (4.47) 



\2L'J (l + 2y)2 \2{m + l)J' 
since the spacing of U' is 

h'^J^^l±^. 

m + 1 m + 1 

The distance y is defined in (|4.34p - (|4.35p . Letting /i — > we have j/ — > (not monotonically), and (|4.47p implies that 



Mi,£,,-4(m+l)^sin^(^^^-^j. (4.48) 
When eN, the end-points of U' arc Xk = —1/2 and Xk+m+i ~ 1/2, and y = 0. o 

Remark IV. 4 When we approximate the continuum by considering e > sufficiently large, we observe that the 
principal eigenvalue fJ-i^u' has the expression (jl.lSp for L ~ 1, in terms of the number ?7i + 2 of oscillators occupying 
the interval U' . Clearly, since < 1 for m > 1 



and we have the bounds 



4 < /ii.c/' - 4(m + 1)-' sin' [ " ) < t:^. (4.49) 

\2(m+ 1)/ 



< ^ — — < Vu < V. (4.50) 



2a + f3 2a + 13 

Besides, for large e > 0, m > 1 is large enough and (|4.49p and (|4.50p justify the approximation 



Ml, (7' ^ 4e sin' 

and the estimation of the power as 



■ 2f ^ \ 
7= , 



4: + Q Hi n> + n 9 / TT \ , , 



^. Numerical results: cubic nonlinearity a = 1 



We now turn to the presentation of the numerical results which starts with the case e = 1. The "cut-off" approxi- 
mation of Proposition IIV. 1 1 takes place on the interval U' of length L' ^ 2 {y = 0.5) and the unit interval U contains 
only one site (to = 1). In Fig. [DJa), the real power of a breather family is plotted using dots against the nonlinear 
parameter a. The lower bound obtained with the "cut-ofF' procedure (j4.45p is shown with a triangle (grey) line. 
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a x„ 

Figure 2: (a) Power of breathers versus nonlinear parameter a in the HDNLS system with cubic nonhnearity [a = 1) and 
e = 1. Symbols (dots) correspond to numerical calculations while the triangles (grey) line represents estimation (|4.45p . The 
continuous (blue) curve corresponds to the estimate (|4.5H) . Other parameters: /3 = = 1. (b) Breather profile (continuous 
(red) curve) against the eigenvector (dashed (black) curve) of (|4.42p - (|4.43p on the interval U' of length L' — 2. The eigenvector 
of (|1.13p in the length L of the chain is represented by the dashed-boxes (blue) curve. 



Notice that it is always below the real power. The qualitative prediction of the pattern of the numerical power as 
given by the theoretical estimate should be remarked, due to the effective approximation of the contribution of the 
linear and the nonlinear part to the power. 

The continuous approximation (j4.5ip in the unit length, plotted with a continuous blue curve, is not satisfied as 
a lower bound for all the values of the parameter a as expected, since we are fairly far from the continuum limit. 
Remarkably, however, we observe that for a quite large regime of the parameter a, the corresponding prediction is 
below the numerical power. This is due to the fact that e = 1 is a critical value for our approximation in the sense 
that for e = 1 the eigenvalue fii^u' in (|4.5ip attains its minimum fii^u' = 4. 

In Fig. [5Jb) the breather profile (continuous (red) curve) is plotted against the eigenvector on U' and the eigenvector 
on the length of the chain L for e = = 1 and /? = 2. Notice that the eigenvector in the length of the system L 
is spread out along the chain (on the scale of the figure it is almost a horizontal line) and its contribution to the 
estimates would be negligible. 

In Fig. Owe present the results of the study for e = 2. Triangles (grey curve) correspond again to the estimate 
(I4.45p . still valid in the interval U' having now length L' ^ 1.414 and the unit interval U contains one site (m = 1). 
We observe the increased quantitative accuracy of the prediction of the actual power (symbols (dots) ) . The continuous 
approximation (j4.5ip in the unit length represented by the dash-dotted (blue) curve is not satisfied as a lower bound 
as predicted by Propositions IIV. l] and ITvT3l Nevertheless, it is worth observing that the continuous approximation is 
only slightly above the actual value. This is connected to the fact that increasing values of e correspond to a closer 
approximation of the continuous limit. The dotted (green) curve below the triangles represents the initial estimate 
(|4.29p with the eigenvalue fii corresponding to the eigenvector of (jl.lSp over the original length L of the system. In 
this case, the estimation of the contribution of the linear part to the power is negligible as (jl.20p shows, thus (j4.29p 
is well below the actual power. 

The effectiveness of the "cut-off" approximation of Proposition lIV.ll and Remark IIV. 21 on length L' , if compared 
with the initial estimate (|4.29p on the length of the system L is even more transparent in the study for e = 3, where 
the results are presented in Fig. |31 In this case L' ~ 1.154 and still m = 1. The curves are traced as in Fig. [31 except 
the new continuous (red) curve which is above the theoretical estimate (|4.29p . This curve corresponds to the lower 
bound in the left-hand side of (|4.5ip . We observe that the prediction of (|4.45p is of excellent accuracy throughout 
the continuation over the nonlinear parameter a and of very good accuracy even versus the nonlinear parameter /3, 
being saturated for large values of /3. It seems that the theoretical estimates capture better the variation over the 
nonlinear coupling coefficient a rather than the onsite nonhnearity coefficient /3. This is due to the fact that through 
the estimation process of Theorem III. 31 the contribution of the hopping nonhnearity is "doubled" by the nonlinear 
coupling with the adjacent sites (see the inequality (|2.22p ). although both nonlinearities are of cubic order in the 
case a = 1. For large values of /3 the manifestation of the power nonhnearity is stronger. More precisely, observe 
in FigHJb) that the convergence of (|4.45p to the real power starts after (3 > 2, i.e. after "doubling" the strength of 
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Figure 3: (a) Power of breathers versus nonlinear parameter a in the HDNLS system with cubic nonhnearity (cr = 1) and e = 2. 
Symbols (dots) correspond to numerical calculations while the triangles (grey) line represents estimation (|4.45|l obtained with 
the the "cut-off" approximation of Proposition |IVjT] The dash-dotted (blue) curve corresponds to the continuous approximation 
(|4.5ip . The first dotted (green) curve from below represents the initial estimate (|4.29|) with the eigenvalue /ii calculated over 
the length L of the system. Other parameters: /3 = f2 = 1. (b) The power and its estimates versus the nonlinear parameter /3. 
Other parameters are chosen as a = Q. = 1. 
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Figure 4: (a) Power of breathers versus the nonlinear parameter a in the DNLS system with cubic nonlinearity (a = 1) 
and 6 = 3. Symbols (dots) correspond to numerical calculations while the triangles (grey) line represents estimation (|4.45p 
obtained with the "cut-off" approximation of Proposition IIV.II The dash-dotted (blue) curve corresponds to the continuous 
approximation (|4.51[) . The first dotted (green) curve from below represents the initial estimate (|4.29|l with the eigenvalue 
calculated over the length L of the system. Other parameters are chosen as /3 = = 1. (b) The power and its estimates versus 
the nonlinear parameter /3. Other parameters are a = Q. = 1. 



the onsite nonlinearity. In this case, the continuous approximation over the unit length approaches further the actual 
power (still, however, from above). 

The approximation procedure considers the cases e = 1,2,3,4, as weak coupling cases, in the sense that the unit 
length U contains only one point and the spacing is ft, > 0.5. Note that U' has different length [L' = 2 for e = 1, 
L' ^ 1.414 for e ^ 2, L' ^ 1.154 for e = 3 and L' = 1 for e = 4.) Since m ^ 1, the eigenvalue in U' given in (|4.44p 
is always ^ii^u' = 2e. Thus, in the weak coupling case, the continuous approximation in the unit length (j4.5ip is not 
valid and (|4.51[) is not satisfied as a lower bound for the power. On the other hand, the discrete approximation with 
the cutoff procedure within (|4.45p becomes progressively better as e is increased. 

Propositions IIV.II and IIV.3I predict that the position of the curves (|4.45p and (|4.5ip should be interchanged when 
e > 4 (ft < 0.5). In this case, the unit interval U contains more than one site (m > 1) and (|4.45p is not vahd. In Fig. 
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Figure 5: (a) Power of breathers versus nonlinear parameter a in the HDNLS system with cubic nonhnearity (cr = 1) and e = 10. 
Symbols (dots) correspond to numerical calculations while the triangles (grey) line represents estimation (|4.45|l obtained with 
the the "cut-off" approximation of Proposition lIV.l] The dotted-dashed (blue) curve corresponds to the estimate (|4.5H) . The 
first dotted (green) curve from below represents the initial estimate (|4.29|l with the eigenvalue calculated in the length L 
of the system and the continuous (red) curve above stands for the estimate (|4.51|l with the lower bound 4e < fj,i{e). Other 
parameters are chosen as /3 = = 1. (b) The power and its estimates are shown versus the nonlinear parameter /3. Other 
parameters a — Q = 1. 



[5]we present the numerical study for e = 10. Here U' has length L' ^ 1.264, the unit interval U contains three sites 
(to = 3) and h 0.316, which can be considered as approaching the continuous limit. Note that for e > 4 we have 
L' > 1, however y — >■ as e is increased. We observe that (|4.45|) is well above the actual breather power in this case, 
while now (j4.5ip provides an adequate approximation especially versus the hopping parameter a. 

Figure [2][b) and Figures [6la)-(b) are showing the breather profiles versus the eigenvectors on U' and the length 
L of the system, and demonstrate the main features of the approximation procedure. A first important feature is 
that both the real breather and the approximating eigenvector for the linear part contribution on U' are localized. 
This is in contrast to the eigenvector associated with /xi (of the original problem) which is extended over the entire 
length L of the system. This approximation of the linear part is effective for values of the weak coupling, where the 
eigenvector on U' has width comparable with the localization length of the breather. In the anticontinuous limit, we 
expect strong localization effects while the eigenvalue fii^u' — 2e becomes negligible again, and the estimates are less 
effective. 

The second feature is that although we are calculating only the contribution to the energy of the sites included 
in U' , the approximation is focusing on these sites being the principal excited ones. Furthermore, Figures [^b) and 
Figures [S] (a), (b) demonstrate a concentration of the "missing" power of the sites outside U' to the most excited sites 
within the eigenvector on U' . This is observable by a comparison of the profiles for e = 2,3,10. From the strong 
coupling to the anticontinuous limit, the breather profile approaches the continuous one, while the eigenvector in 
U' converges to the continuous eigenfunction. Then, both the difference between the breather and the eigenvector 
width as well the difference of their "peaks" becomes constant, and again, the estimates are becoming less effective. 
Besides, the methods of this paper are making use of the properties of the discrete phase space and should be 
extended appropriately in function spaces in order to capture effectively the behavior of the continuous counterpart. 
Nevertheless, in this setting of larger e, the continuum variant of the approximation over the interval U' yields a 
suitable lower threshold for the breather power. 

3. Numerical results: Quintic nonlinearity a = 2 
Concerning the case of non-cubic nonlinearity (ct 7^ 1), an explicit estimate from equation (|4.28p comes out 

<V, o->l, A^>1, (4.52) 




{2aN)^ (T- 1 
(/3a) ^ ^ , 
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Figure 6: (a) Breather profile for e = 3 (continuous (red) curve) against tfie eigenvector (daslied (black) curve) of (|4.42p - (|4.43p 
on the interval U' of length L' ~ 1.154. The eigenvector of (|1.13|l in the length L of the chain is represented by the dashed-boxes 
(blue) curve. Other parameters are a = 1, /3 = 5, = 1. (b) Breather profiles for e = 10. Here L' ~ 1.264. Other parameters 
a = 1, /3 = 5, f7 = 1. 



with some restriction on the parameters a, Q, and e given in Theorem III. 31 C: 
(i) for all 17 > if 



and (ii) for all e > if 



e > 



VL > 



{2aN)'- 



1 



Aicr 



{2aN)^ a - 1 



cr > 1, > 1, 



cr > 1, iV > 1. 



(4.53) 



(4.54) 



(/3a) — 

In the non-cubic case, the cut-off approximation of Proposition IIV. 1 1 leads to 

Corollary IV. 1 Let a > \, N = \. Then the estimate 114 -5 2^ is valid with fj.i replaced by 

A. fii^u' = 4esin^ { 2(^+1) ) °-^y e > 0- 

B. fjLi^u' ^ 4esin^ (^Te) e > is sufficiently large. 

The theoretical estimates proposed in section IIIII for infinite lattices can also be used. While an infinite lattice 
cannot be modelled numerically, the estimates of section IIIII can serve as alternatives to those summarized above for 
the finite lattice. The unspecified parameter z^crit involved in p.2ip . in the estimate p.23p and restrictions p.24p - 
p.25p has been determined by justified heuristic (and rigorous in the case of "large" cr) arguments in [111, Section III, 
pg. 7]. For instance it was revealed that the value i^crit = 1 is valid for all > 1 and cr > 1. Furthermore, this value 
is of very good accuracy for N = 2 and excellent for iV = 3. Let us also recall that this value covers when cr e N, the 
cases which arc of main physical interest (see also [s^ considering integer values of cr > 2/N). 

For j/crit — 1, Theorem IIII.2I predicts that for supercritical nonlinearity a >2lN any breather solution must have 
power 

72 



-^crit < 



^crit is the positive root of the equation 

2aNR'^ 

Theorem IIII.4I gives explicitly 



PR 



2<j 



AeN 
2a + I 




{2aN) 



<T-1 (J 



1 



2ct-^ 1 



(Pa)-' 



< R^, cr > 2 when iV = 1 and cr > 1 when iV > 2, 



(4.55) 



(4.56) 



(4.57) 
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in either of the cases below: 

(i) for all > and lattice spacing satisfying 



e > l)(2a + l) ^ ^ ^ ^^^^ iV = 1 and a > 1 when iV > 2, (4.58) 

(/3a)^ 4iVa ' - 

and 

(ii) for all e > and frequencies 

{2aN)^ cr - 1 

Q. > 4 , cr > 2 when N =1 and cr > 1 when TV > 2. (4.59) 

i/3a)— ^ 

Additionally other choices of the parameter e in the Young's inequality trick (see Theorem IIL3C ). give versions 
of the estimates valid with different restrictions on the coupling parameter e or the frequency fl. Together with the 
choice used in Theorem III. 3C another interesting one is the standard e = 1 corresponding to the version of p.23p 



^ ^ AeN (cr - l)(2a7V)^ 



(7/3+1 \ 2(7+1 



The estimate (|4.60p is valid 

(i) for all 17 > and lattice spacing satisfying 



< R^, a >2 when iV = 1 and cr > 1 when iV > 2 (4.60) 



(2aN)~(a - l](2a + 1] 
e > ^—^ ^ '-^ cr > 2 when iV = 1 and a > 1 when iV > 2, (4.61) 

47Vcr ' - - ' ^ ' 



and in the case 

(ii) for all e > and frequencies 



(2aN)~{(7 - 1) 

n > ^ '- a>2 when TV = 1 and (t > 1 when N>2. (4.62) 

cr 
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Figure 7: (a) Power of breathers versus parameter /3 for supercritical nonlinearity cr = 2 in tire HDNLS system. Symbols 
correspond to numerical calculations, the dash-dotted (blue) line represents estimation (|4.57p and the dashed line (green) 
estimation (|4.60|l . The estimate (|4.52|I -A. of Corollary II V . II corresponds to the triangles (grey line), and the estimate (|4.52I) -B. 
with the continuous (red) line. Parameters: a = 0.01, SI = e = 1. (b) Power versus frequency Q for supercritical nonlinearity 
cr = 2. Parameters: a = 0.5, 13 — 5 and e = 1. 



Regarding the quintic nonlinearity, more specifically, we have performed a test of the estimates (|4.52p . (|4.57p and 
(|4.60p by fixing cr = 2, and e = 1. With these choices, restrictions (|4.54p - (|4.59p and (|4.62p reduce to the very simple 
conditions il > 0^/(3 and fl > 2(2? . We expect all the estimates to be satisfied as thresholds due to the increased 
strength of the power nonlinearity absorbing the contribution of the linear part, even in the case of (I4.52l) -B.. which 
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is not justified theoretically. In FiglTl^a) we have plotted the estimate (|4.52[) -A. of Corollarv lIV.il with triangles (grey 
curve), and its case B. with the continuous (red) curve. The dash-dotted (blue) and dashed (green) lines correspond to 
(|4.57|) and (|4.60|) respectively. The numerical power (symbols) was obtained varying /3 for a small hopping parameter 
a ~ 0.01 and = 1. Note that all the estimates are good, although (|4.60p is better than (|4.57p for large /3 while 
(|4.57p behaves better when /3 < l/cr. 

In Fig. Ill^b) we have plotted the breather power against f2 choosing /3 = 5 and a = 0.5. Condition ()4.54|) is fulfilled 
for f2 > 0.05 and condition (|4.62p is fulfilled for J7 > 0.5. In the latter region, since /3 is quite large, the estimate 
(|4.60|) behaves clearly better than (|4.57p . It is interesting to realize that (|4.52p is worse than (|4.60p for large enough 
frequencies. 



B. Defocusing case (a < 0, /3 =< 0) with Dirichlet boundary conditions. Solutions tp„{t) = e '^'^n, H > 0. 



In the defocusing case the results on the theoretical estimates are restricted to frequencies ft > ANe. In this case, 

setting for convenience k ~ —a > 0, A = —(3 > 0, the results of Theorem III. 41 state that for all cr > the lower bound 
for the power of the staggered breathers is given by the positive root R^,^d of the equation 

Ax^" + 2kA^X^ - (O - 4e7V) = 0, a > 0, N > I, r2 > 4eiV, (4.63) 

and the power of staggered breathers satisfies 

Rl a < r, for all CT > 0, iV > 1, n> AeN. (4.64) 

In the defocusing case and cubic nonlinearity, the lower bound for the power is 



2kN + X 

The explicit estimate valid for ct > 1 is 



<r, n> AeN, N >1 (7 = 1. (4.65) 



2A \^ (Aa)— 



(2kN)^ cr — 1 

<V, cr>l, N>1, n>4eN+^ . (4.66) 

(Act)— cr 



The results of the numerical tests in the defocusing case are similar to those of the focusing case and can be 
summarized in the following points: 

• The theoretical estimates are always below the numerical power and approximate quite well the nonlinear part 
of the contribution to the power. 

• The lower bound (|4.64p is always above the explicit estimate (|4.66p 

• Estimate (|4.66p behaves better for small values of the hopping parameter a and large exponents a. 

These observations are corroborated by the results of FiglSl Squares and the upper continuous curve correspond 
respectively to the numerical power and estimate (|4.64l) for /3 = —1 and a hopping parameter a = —0.5. The estimate 
becomes much closer the real power fixing /3 = —5 and a = —0.01 (see pluses and the lower continuous curve). 



V. CONCLUSIONS 



In the present work, we generalized the considerations of energy thresholds in the setting of a DNLS model with 
generalized nonlinear (Hamiltonian) hopping terms. Different types of bounds were provided for the power both for 
finite and for infinite lattices, by using appropriate estimates for the linear coupling and nonlinear hopping terms. 
A fixed point method establishing the contractivity of an appropriately defined operator was also used to establish 
that for a given parameter set, there is a critical power, below which it is not possible to sustain such nonlinear 
waveforms. Finally, some dimension-dependent estimates were given based on the interpolation inequality of the 
Gagliardo-Nirenberg type, in a spirit similar to the work of [2^ . 

Further improvements of the main theory have been considered and proved, appreciating the interplay of the 
nonlinear and linear term contributions within the true solitary wave solutions, taking into account their spatial 
localization. The obtained bounds were tested numerically and in all the cases where the theory was expected 
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Figure 8: Power of breathers versus parameter a in the defocusing case. Squares (pluses) correspond to the numerical power 
found for /? = — 1,q = —0.5 (/3 = —5, a = —0.01) while the upper (lower) continuous line represents estimate (|4.66|l . Other 
parameters: Q. = —2, t — 0.25. 



to be applicable, it was found that the numerical solutions satisfy the predicted norm inequalities. This aspect also 
provides details on the parameter regimes (weak linear coupling) which tend to saturate the corresponding theoretically 
obtained bounds. 

We are leaving as an interesting open direction for a future work, to examine the behavior of the energy bounds 
when the size of the lattice is varied. This question is taking into account the effect of the transition from finite to 
infinite lattices, on the localization properties of the solutions. This task could be based on a generalization and use 
of the machinery developed in [32j . as well as, of the relevant localization estimates. Such a generalization could be 
of particular interest, in the case of multidimensional lattices. 

It would be also interesting and relevant to examine how corresponding bounds can be generalized to other classes 
of models, including ones of the nonlinear Klein-Gordon or FPU type (or mixed ones), incorporating different types 
of onsite and intersite nonlinearities. Especially useful, albeit arguably more difficult, to extend the main strategy to 
continuous models. Such tasks will be considered in future publications. 
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